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We study a mixture of one-dimensional bosons and spinless fermions at incommensurate fill- 
ing using phenomenological bosonization and Green's functions techniques. We derive the relation 
between the parameters of the microscopic Hamiltonian and macroscopic observables. Galilean 
invariance results in extra constraints for the current current interactions. We obtain the exact 
exponents for the various response functions, and show that superfluid fluctuations are enhanced by 
the effective boson-fermion density-density interaction and suppressed by the effective boson-fermion 
current-current interaction. In the case of a bosonized model with purely density-density interac- 
tion, when the effective boson-fermion density-density interaction is weak enough, the superfluid 
exponent of the bosons has a non-monotonous variation with the ratio of the fermion velocity to the 
boson velocity. By contrast, density- wave exponent and the exponent for fermionic superfluidity are 
monotonous functions of the velocity ratio. 



I. INTRODUCTION 

The recent development of atom trapping technologies 
0, has led to the opening of new research fields of 
many-body physics. In particular, the possibility of con- 
trolling to a large degree the shape of the trapping poten- 
tial has permitted the experimental realization of nearly 
isolated low dimensional interacting boson systems jj, 
and the observation of the bosonic Mott insulator jj. 
Moreover, the use of Fano-Fcshbach resonances [B|, @| al- 
lows experimentalist to control interaction strength in 
the low dimensional atomic gases 041] • This has per- 
mitted the observation of the Girardeau [To| fermion- 
ization of one-dimensional bosons with strongly repul- 
sive interactions [ll], EH ■ From the theoretical point of 
view, fermions as well as bosons with repulsive interac- 
tions in one dimension are in the Tomonaga-Luttinger 
liquid state (l3l - [l6j . The bosonization technique [l]| al- 
lows the description of their low-energy physics in terms 
of collective phonon-like modes as well as the computa- 
tion of the various correlation functions. It is found that 
in the Tomonaga-Luttinger liquid state at zero tempera- 
ture, only quasi-long range order (whether superfluid or 
density wave) can be obtained, with correlation functions 
decaying as power laws, the exponent of the power law 
being a rational function of the Tomonaga-Luttinger pa- 
rameter K. This parameter K depends on interaction. 
In the case of bosons with contact interaction [l?}, this 
parameter can be calculated exactly, and it is found that 
K > 1 making superfluidity the dominant correlation. 
When the interaction becomes extremely repulsive, or 
the density becomes very low, the Girardeau model is re- 
covered [To] and the exponent K goes to 1. In the case of 
spinless fermions, perturbative calculations (l8l - l20j show 
that the exponent K is smaller than 1 for repulsive in- 
teractions (thus favoring density waves) and larger than 
1 for attractive interactions (thus favoring a p-wave type 



superfluidity). For finite temperatures, the correlation 
functions decay exponentially, with a correlation length 
inversely proportional to temperature. Such behavior 
gives rise to susceptibilities that diverge as a power law 
with temperature, with an exponent determined by the 
zero temperature correlation exponent. Recent experi- 
ments [2ll425l | with trapped atoms have shown the pos- 
sibility of realizing many-body systems mixing fermions 
and bosons. From the theoretical point of view, mix- 
tures of bosons and spinless fermions in one dimension 
were considered in Refs. [2fj| - [28l It was found that suffi- 
ciently strong repulsion or attraction between fermionic 
and bosonic atoms could result in respectively a phase 
separation or a collapse, while interactions of moderate 
strength would lead (in the case of incommensuration 
between the atom densities) to the formation of a two- 
component Tomonaga-Luttinger liquid. In Refs. (2?1 HH it 
was proposed to describe the two-component Tomonaga- 
Luttinger liquid in terms of polarons, and various ex- 
ponents of correlation functions were obtained. An in- 
tegrable model of interacting bosons and fermions was 
considered in Refs. [29l - [33l in which a two-component 
Tomonaga-Luttinger liquid was found. The mixture has 
also been investigated in numerical calculations [33 - 140} . 
More recent experimental work has considered mixtures 
of fermionic atoms. Motivated by these experiments, 
there has been some numerical studies of mixtures of 
three fermionic species fill ]. 

Although the one-dimensional boson fermion mixture 
has been studied extensively including analytical expres- 
sion of the exponents [13, [HI of the equal time correlation 
function, the previous work has emphasized the polaron 
correlations over the density wave correlations. It is thus 
worthwhile to analyze the competition between density- 
wave and superfluid correlation as interactions in the sys- 
tem are varied. Besides the question of the effect of in- 
teraction on the exponents, an important question rele- 
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vant for experiment is the evolution of the exponents as a 
function of relative concentration of fermions and bosons 
in the mixture. Fermions could have two opposite ef- 
fects. On the one hand, fermionic particles can give rise 
to retarded attractive interactions between the bosons, 
which can compensate the direct boson-boson repulsion 
and favor superfluidity and possibly collapse of the mix- 
ture. On the other hand, fermionic particles do not form 
Bose condensates, and may thus disfavor the formation 
of a superfluid state. For the simple Bose-Fermi mixture 
model, it is useful to examine the effect of the doping on 
the correlated states of respective particles, e.g. , how the 
superfluidity is enhanced or suppressed by the interaction 
between boson and fermion. A last question, in order to 
compare numerical and analytical result, is to develop a 
phenomenological bosonization scheme allowing for the 
calculation of the parameters of the low energy effective 
Hamiltonian from ground state energy calculations. 

In the present article, we present a study of mixtures 
of fermionic or bosonic atoms using bosonization and a 
Green's function formalism. After discussing briefly in 
Sec. |n] the phenomenological bosonized Hamiltonian of 
the model and its relation with the macroscopic observ- 
ables, we focus in Sec. IIIII on the calculation by Green's 
function techniques of the exponents for superfluidity and 
density waves. We find a relation between the products 
of density wave and superfluid exponents for the bosons 
and the fermions. We also give an expression of the cor- 
relation lengths at finite temperature. In Sec. IIV1 we use 
the Green's function technique to obtain the Bragg scat- 
tering intensity following Ref. |42|. Then, in Sec. [V] wc 
consider the effect of the variation of the density of one 
of the two atom species first in the case of the model of 
Ref. H^. We show that for this model, the variation of 
the bosonic superfluid exponent with fermion density is 
not monotonous in the case of a weak interaction. Wc 
also consider the effect of the effective current-current 
interaction and show that it suppresses superfluidity. 



When the boson and fermion densities are not com- 
mensurate with each other, the bosonized Hamiltonian 
reads MfM: 



H 



E 

a , b 



dx 
2^ 



[7r 2 A/ a6 n a n fc + N ab d x (f> a d x ^ b ] 



(2) 



where [(j> a (x),Hb(x')] = iS a ,bS{x — x') and a,b € {B, F}. 
Indeed, the Hamiltonian ((TJ) is invariant under a parity 
transformation. Under such a transformation, 4> a {x) — > 
—4>a{—x) and TL a (x) — > — n a (— x), so that quadratic 
terms H a d x ipb change sign under parity, and Eq. J2]) 
is the most general Hamiltonian quadratic in n a and 
d x (j> a for a = B,F. The matrices M a b and N a b in Eq. 

are real symmetric and are deduced from variations 
of the ground state energy Eos of a finite system of 
size L with (respectively) changes of boundary condi- 
tions ip a {L) = e l(Pa tp a {0) and changes of particle densities 
Pa = N a /L: 



M„ 



N ab = 



d 2 E GS 

dlfiadlfb ' 

1 d 2 E GS 
it L dpadp b 



(3) 
(4) 



In the case of a Galilean invariant Hamiltonian such as 
Eq. (JTJ) , it is possible to further constrain the elements of 
the matrix M with the relations (see Add. |A"|) : 



irp F = MffM f + MbfMb, 
irp B = MbfMf+MbbMb- 



(5) 
(0) 



So that knowing just one of the parameters Mff, Mbb, 
and Mbf fully determines the matrix M. If one starts 
from the limit Vbf = and applies bosonization first, 
and then reintroduces the interaction Vbf(x — x') = 
VbfS(x — x 1 ) as a perturbation, to lowest order one ob- 
tains the Hamiltonian: 



II. HAMILTONIAN 

We consider a mixture of spinless fermions and bosons. 
The Hamiltonian reads: 



H = 



dx 



-^B d l^l 



i 

2 

Vbb(x 



2Mb 2M f 
dxdx'[VFF(x — x')pf{x)pf{x') 
x')pb{x)pb(x') + 2Vbf(x 



E 

v=F,B 

Vbf 

IT 2 

in which u u K u 



dx 
2^ 



11 

u v K v (ttU„) 2 + -^-{d x <t> v ) 



dxd x (j) B d x 4> I 



(7) 



.,10) 



as a result of Galilean invari- 



ance of Eq. (fTl) , -fe- = — tA— — where pi^ is the density 

of particles, and k„ is the compressibility. In the case 
of non-interacting fermions Vff = 0, one has Kp = 1. 
x')pB(x)pF(x(J\f 0Y hard core bosons [Io|, Kb = 1, while for the Lieb- 
Liniger model [13], Vbb(x — x') = V B bS{x — x'). one has 
Kb > 1, with K B — > oo when V BB — > and K B — > 1 
for Vbb — > oo- Note that at this low order in pertur- 
bation theory, no term Mbf is present. This term is 



where Mf,Mb are the masses of (respectively) the 
fermionic and bosonic atoms, Vff,Vbb,Vbf (respec- 
tively) the fermion-fermion, boson-boson and fermion- 
boson interactions, and ipF,ipB the (respectively) bosonic 
and fermionic annihilation operators. We have also de- 
fined the density operators: p v = The model ([JJ 
can be treated by the method of bosonization [l4| ■ 



expected to appear in second order perturbation theory, 
along with corrections to the bare terms Mff and Mg b ■ 
A Hamiltonian with quadratic interactions similar to Eq. 
@ but also comprising interband tunneling terms was 
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considered in the context of a two band model of inter- 
acting spinless fermions 43]. A Hamiltonian equivalent 
to Eq. ([7]) has been studied by path integral methods 
as a model of one-dimensional electrons interacting with 
acoustic phonons [3,|45j]. Due to the quadratic character 
of the Hamiltonian ([2]), its spectrum is readily obtained 
[26j as two branches uj±(q) = u±\q\ with linear spec- 
trum, showing that its ground state is a two component 
Tomonaga-Luttingcr liquid provided the velocities u± of 
both components are real. The vanishing of the veloc- 
ity u_ is an indication of instability [2t| either towards 
phase separation (in the case of repulsive boson-fermion 
interaction) or collapse (in the case of attractive boson- 
fermion interaction) . 



III. SUPERFLUID AND DENSITY WAVE 
CORRELATIONS 

From the diagonalization of the Hamiltonian ©, 
bosonization allows to obtain the exponents of the vari- 
ous correlation functions [28j . In this paper, we will use 
a different (but equivalent) approach to compute the cor- 
relation functions of the system. 

By an equation of motion method, we first obtain the 
Green's functions: 

G ab (x,r) = -{T T (Mz,T)-M0,0))<h(0,0)), ( 8 ) 
G ab (x, t) = -(T T (6 a (x, t) - O (O, 0))0 6 (0, 0)), (9) 

where 6 a = n J dx'U a (x') and a, b can be B or F . From 
these Green's functions © and ©, we find the correla- 
tion functions of exponential fields as: 

T^E^^^e^E^^ ' ^ = e^^" XaXbGab{ \lh) 

where Xb and are real numbers. Such method is of 
course applicable to cases with more than two compo- 
nents, as long as the Hamiltonian remains quadratic in 
the fields (j) a and 8 a . 

To derive the Green's functions, we start from the 
equations of motions in Matsubara time of the fields <p a 
and n a read: 

d T 4> a (x,T) = [H,(t> a ]=-mY^ M ab*lb(x,T), (12) 

b 

9 r n Q (.x,r) = [H,U a ] = --Y J N ah dl<j> b . (13) 

TV ■ * 



Going to Fourier space, we obtain (a, b = B,F): 

G ab (q,u) n ) = -^((u 2 n + (MN)q 2 )- 1 M) ab , (16) 
where lo„ = 2nnT. From Eq. (|TB)) : 

(17) 

where we have introduced the cutoff a. By using the 
duality transformation d x (j) a = ^P a , d x 9 a = Ttll ai we 
obtain equations of motion for G ab similar to Eq. (|14[) 
with the roles of M and N exchanged. Thus, in Fourier 
space, we have: 

G ab (q,u n ) = -ir((u% + (NM)q 2 )- 1 N) ab . (18) 

The expressions (fll)|) and (fl"8|) show that the retarded 
Green's functions have poles for iui equal to u±\q\ where 
u 2 - < u 2 ^ are the two eigenvalues of MN. Stability re- 
quires that u± > and thus det(A/7V) > 0. Since we 
know from App. [5] that dct(M) > in a Galilean in- 
variant model, instabilities occur for det(Af) = 0. From 
the definition of N, Eq. ([1]), such instabilities are either 
collapse or phase separation. 

Further, for zero temperature and r = 0, we can obtain 
a general form for G using Eq. (fT7|) . We find: 



G(x,0) = -{MN)-V 2 Mhx 



G(x,0) = -(NM)- 1/2 N\n 



\J x 2 + a 2 



V 1 x 2 + a 2 



a 



Then equations (fTO]) and (fTTj) lead to: 



\Jx 2 + 1 



\/x 2 +, 



(19) 
(20) 

^A(J\/JV) _1/2 A/A 

( 



(21 



where *A = (Xb,Xf)- We can define the matrices 77^ = 
(MN)-^ 2 M and r)e = (NM)~ X / 2 N = (NMf^M' 1 . 
These matrices yield the exponents for the exponential 
fields. One can see that the duality relations become 
V^Ve = Ve tr lct> = 1- For the two component system, one 
has the identity: 



The equations of motions for the Green's functions 
G ab (x,r) thus read: 

d T G ab (x,r) = «v][]M oc (T T n c (z,T)^(0,0)}, (14) 

c 

d 2 T G ab (x,T) = n5(x)S(T)M ab -J2( MN )acdlG cb (x,T). 

(15) 



(MN) 



-1/2 _ 



■(MN)- 1 , (23) 



which can be checked by applying the right hand side of 
the formula to each eigenvector of MN. We thus have: 



1)0 



(M + U+U-N- 1 ) 



(N + u+U-M' 1 ) 

u + + u_ 



(24) 
(25) 
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In the case where det(iV) — > 0, we see that ij^ will 
have matrix elements going to infinity as uZ , whereas 
r/e — > N/u + . Therefore, near a collapse or a phase sep- 
aration the density wave exponents are divergent, while 
the supcrfluid exponents remain finite. 

For nonzero temperature, the sum Eq. (|17[) is domi- 
nated for long distances by the term with n = 0. One 
finds that: 

\ |G(^o=0)r-l)^ = \t\x\N~\ (26) 
i J G( g ,u, o = 0)(e^-l)g = \t\x\M~\ (27) 



where we have defined pf,o( x ) = (ifj + Y (x)ip + (x) + 
{ip~)' { {x)ip~ (x). For the bosonic density wave fluctua- 
tions, using the Haldane expansion of the density, 

p B {x) = p$ - -^d x (f>B + ^2 A m cas(2m4>B ~ 2imr p^ x\34) 



we find that the density density correlation function of 
the bosons reads: 



{p B {x) PB {Q))^{pff 



(B) 

Vadw 



E 

m>l 



cos(2rmrpgx) 
Att 2 x 2 ^ 4^ (xlaY^^Aow 



(35) 



so that the correlation functions decay exponentially for 
long distances, 



J2 a K<Pa(x,0) -iJ2 a A °<M0,0)\ 



■W,(28) 



,(29) 



The dominant correlations are at wavevector 2itp^) and 
are characterized by the exponent: 



(B) _ 

^Iadw — 



M 



BB 



U+U— 

dct(AT) 



N, 



FF 



(36) 



with thermal correlation lengths given respectively by 
&(A) = 2/(^T*A(iV)- 1 A) and &(A) = 2/(7rT*A(M)- 1 A). 
We note that near the instability, the correlation length 
£</>(A) goes to zero, in accordance with the reduction of 
density wave ordering found at zero temperature, while 
the length £g remains finite. 



A. Atomic density wave correlations 

To characterize the atomic density-wave (ADW) order- 
ing of the fermions, the field operator of which is given 
by: 



i) F (x) 



1% 



V27ra 



(30) 

where k F — irp)^ ' is the Fermi wavevector, a is a short- 
distance cutoff, and ± label the two Fermi points, we 
have to calculate the correlation function of the 2k F 
component of the density operator if; 1 ' (x)ip(x), namely 



(P2k F ~ {tp^ip F ): 
(pL^)^ F (0))~(e-* W+l2fe(0 ») 



( F) 

The density wave exponent fl^nw 
2 



(x/a) 11 AD 



w(31) 



(F) 
'I ADW 



U+ 



M 



FF 



U+U- 

det(N) 



N 



BB 



(32) 



For Vbf — 0, the exponent of Eq. (I3"2l reduces to 2Kp. 
Near the collapse (for Vbf < 0) or the phase separation 
(for V B f > 0) which is obtained at it_ — > 0, we note 

that the exponent ^Iadw IS diverging as ~ l/it_. The 
fermionic density-density correlation at small wavevec- 
tors can also be obtained from bosonization. We have: 



(PF, (Z)PF,0(0)) = {pff 



(F) 
47T 2 X 2 ' 



(33) 



This exponent is obtained from the fermionic exponent 
by the exchange (M B b,N ff ) «-> (M ff ,N B b)- The ex- 
ponent is also divergent when u_ — ¥ 0. 

It is also interesting to consider the cross correlations 
between bosonic and fermionic density. One has: 

(o) (o) M B f - u+u-N B F/det(N) 1 



( P b(x)pf(Q))=P K b'pF- 



2ir 2 x 



so that the non-uniform components of the densities of 
bosons and fermions remain uncorrelated. The cross den- 
sity correlations vanish when Vbf = and are positive 
when Vbf < as a result of the boson-fermion attrac- 
tion. 

For finite temperature, the correlation functions of 
Eqs. (|3"Tj) and (|31)|) decay exponentially, with correlation 
lengths given by: 



(F) _ 2det(JV) 
Sadw- nNBBT > 

2det(7V) 



AB) 
?ADW 



ttN ff T ' 



(38) 
(39) 



B. superfluid correlations 

Fermions, due to their spinless character, can only 
present p-wave type superfluidity. The order parameter 
is —iip(x)Vi/}(x) and can be expressed using the decom- 
position of Eq. (|3T)]) in the form: 



^f( x )^f ( x ) 



2na 



(40) 



The order parameter for p-wave superfluidity exhibits al- 
gebraic correlations: 

(^MV^XV^W^O)) 1 ) ~ {x/a)-^\ (41) 
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with the exponent for fermion superfluidity: 



(F) 

Vs = 



N. 



FF 



U+U- 

det(M) 



M 



DB 



(42) 



We now turn to the superfluid fluctuations of the bosons. 
The quasi-long range superfluid order is characterized by 
the correlation function: 



(il>Ux)M0)) ~ (e-^W+^W) ~ (x/a)~ 
Our result for the superfluid exponent rj s is: 
1 



(B) 

Vs 



2(u + +u-) 



N, 



BB 



-M. 



FF 



(43) 



(44) 



det(Af) J 

The superfluid and the density wave exponents 
are not independent from each other. Indeed, not- 
ing that ?i + it_ = det(M.ZV) 1 / 2 , we have that: 
4rj { s B) = (det(7V)/det(Af)) 1 /2 7? (f z ) H/ and ^ = 

(det(A^)/det(M)) 1 / 2 

Vadw This i m P nes that the expo- 
nents satisfy the relation VadwVs = ^iadw^S ' 

Turning to the finite temperature case, the thermal 
lengths are: 



2det(M) 
ttM bb T ' 
2det(M) 
irMppT ' 



(45) 
(46) 



IV. BRAGG SCATTERING 

According to Ref. 42j, the Bragg scattering intensity 
is proportional to the imaginary part of the retarded 
density-density response functions x(g,a>). Retarded 
density response functions Xab(Q>w) i a , b = B,F) can 
be obtained from the Fourier transform of Green's func- 
tions HO) as: 



G ab {q,uj n )\ 



iu) n — >uj-\-iOj r ■ 



(47) 



Using the expression of the Fourier transform (|I6[) from 
Sec. IIIII we obtain for cj > : 

U-{M ab -u\{N- l ) ah ) x/ 
2 — °\ u ~ u -\l\) 



ImXab(q,u) = — — 



u+iuliN- 1 )^ - M ab ) 



S(uj - u+\q\) 



(48) 



Equation (|48l) predicts peaks at frequencies u±|<z|- The 
matrices M and N can be deduced from the spectral 
weight of these peaks. One has the following sum rules: 



f°° M 
/ dujlmxab{q,u) = —{rj^at, 

Jo 2 TT 





00 duj 



imxab{q,w) 



-AN- 



jab: 



du)(jImXab(q,u) = —M ab 



(49) 
(50) 
(51) 



The first sum rule is simply a restatement of our deriva- 
tion of the equal time Green's function in Sec. HH The 
second sum rule is a Kramers-Kronig relation giving the 
real part of the zero frequency (matrix) density-density 
response function as an integral of its imaginary part. 
Since the real part of the density density response func- 
tion is the compressibility, the second result is not sur- 
prising. The last sum rule is a consequence of cur- 
rent conservation. Indeed, using current conservation, 
one can relate the density-density response function to 
the current-current response function. Using again a 
Kramcrs-Kronig, the last integral is shown to be equal 
to the static current-current response function. 



V. VARIATION OF EXPONENTS 

In this section, we first consider the model of Eq. ([7]) 
with Mgp = 0. For the two-component case (a, b g 
{B,F}), the velocities are found as [261 ]: 
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± 



v, 



BF 



ubKbufKp[52) 



Using the Green's function methods of Sec. [TT1 we obtain 
the following expressions for the exponents: 



(B) 

Vadw 



(F) 

Vadw 



(B) 

V s 



(F) 

Vs 



ubK b 

U + U- 

ufKp 

I 

Au B K B 
I 

upKp 



,2 1 



+ U- + 



U + + U- 



U + + U- 



(53) 
(54) 
(55) 
(56) 



Finally, using that u 2 . 



2 „,2 



IlpM 



F U B 



(Vbfu) 2 /it 2 where u 



l f + u B and u + u_ = 

= \JupKpubKb-, we 
can express the exponents entirely as functions of 
uf,ub,Kf,Kb and Vbf- Expanding to second order 
in Vbf we find: 



(B) 

Vs 



(F) 

Vs 



(B) 

Vadw 



(F) 

Vadw 



1 



VgpK F U F 



2Kb 
2 

~K~f ~ 



= 2Kb 



2Kt, 



4tt 2 ub(uf + ub) 2 
Vbf k bu b 
tt 2 uf(uf + ub) 2 

V^ f K f Kb{2ub+u f 



+ 0(V£ F ), 



(57) 

+ ()(V£ F ), (58) 

+ o(uj F ) , 



2tt 2 {uf + u b ) 2 ub 

(59) 

" L;) +o(yi F ) ■ 



V 2 f K f K b (2uf 



2ir 2 (up + ub) 2 uf 



(60) 



This shows that superfluidity is enhanced by the boson- 
fermion interaction. The enhancement of fermionic su- 
perfluidity becomes weaker as up is increased. By con- 
trast, when up increased, the enhancement of bosonic 
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FIG. 1: The contour plot of 2K B r) ( f ) with the fixed £ = 
uf/ub = 2 on the plane of N and M(< 0) which are the di- 
mensionless quantities defined by N = Nbf\/KbKf /^/ubUf 
and M = Mbf(VKbKf)^ 1 /^ubuf, respectively. For 
N = (M — 0), the quantity 2Kbvi^ increases (decreases) 
monotonously as \M\ (\N\) increases. 



superfluidity is non monotonous, being maximum for 
up = ug. At the same time, the exponent of density 
wave correlations for the bosons and fermions are in- 
creased by Vbf, indicating a reduction of the density 
wave quasi-long range order. 

We note that the density-wave exponents Eqs. (|53|) 
([54")) can be rewritten respectively as = (u% + 

u+uJ) / (2ubKb\J u B + up + 2u + u-) and r)^ = (up + 
u+U-) /\upKf\/u B + Up + 2u+U-], which are both in- 
creasing functions of Since from Eq. (|5"12|) . is 
a decreasing function of |Vrf|i the superfluid exponents 
are decreasing functions of |Vbf|- In other words, boson- 
fcrmion interactions lead to an enhancement of superfluid 



fluctuations. 

In the context of electrons coupled to acoustic phonons, 
an analogous enhancement of fcrmionic superfluidity was 
already noticed [45] . The interpretation of the result (|42)l 
is that the sound modes of the bosonic fluid interact with 
the fermionic atoms exactly as the acoustic phonons in- 
teract with the electrons in the model of Ref. [45]. The 
interaction of the acoustic modes with the fermions gives 
rise to an effective attractive interaction that enhances 
superfluidity. The enhancement of bosonic superfluidity 
can also be understood in terms of effective attraction be- 
tween bosons generated by the sound modes of the Fermi 
fluid. We note that in the context of electronic systems, 
the collapse or phase separation that results from exces- 
sively strong Vbf interaction is known as the Wentzel- 
Bardeen instability [46], |47| . The enhancement of super- 
fluid correlations in mixtures could be tested by varying 
the strength of the boson- fermion interaction Vb f in a 
40 K- 87 Rb mixture using a Feshbach resonance Su- 
perfluid correlations will present a minimum when Vbf 
is tuned to zero by interaction. 

Next, we discuss the case in the presence of both 
Nbf{— Vbf 1^) and Mbf in Eq. ©, which give the 
Hamiltonian, 



H 



E 

v=F,B 
N B F 
7T 



dx 
2^ 



UvKJ-kT1 u ) + -^-(d x (j> v y 
iW 



(61) 



dxd x 4>Bd x <pF + itMbf / dxHs^F- 



As shown in Appendix [Aj we have the constraint Mbf < 
from the Galilean invariance of the system (see Eq. 
(lAlpp ). so we restrict ourselves to the case of Mbf < 
in the following. The velocity is obtained as 



Up + III 



NbfMbf ± 



ubK b u f Kf N B f 



MbfUb \ 
K B Kpu F ) 



N, 



BF 



MbfUf 

K b K f ub 



(62) 



Thus the excitation is stable for 



From Eq. (|44[) , the exponent for the superfluidity for bo- 
son is calculated as 



uW_ = (u B u F - M 2 bf I{KbKf)){ubu f - N B pK B K F ) > Qf63) 



2K, 



l/(2u B K B ) 



y/u 2 B + up + 2N bfMbf + 2u+m_ 



1 + + 2£NM + 2£J (1 - N 2 )(l - M 2 ) 



u b u f {ubuf - N B pK B K F ) 



1 + 



1-M 2 



(64) 
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where N — Nbf\/ KbKf / \Jubuf and M = 
MbfW KbKf)^ 1 I 'y/UBUF- The second line shows that 
2K B r]g depends on N, M, and £ = u f /ub- The sta- 
ble conditions for u± are expressed as —1 < N,M < 1. 

~ ( b) 

At the boundary of \M\ = 1, rj s becomes infinite while 
it stays constant at \N\ = 1 (\M\ ^ 1). Since TV (M) 

f B) 

has an effect of decreasing (increasing) r/ s , they com- 
pete with each other. In Fig. [TJ the contour plot of the 

( B) ~ 

quantity IKb^s is shown on the plane of N (horizontal 
axis) and M (vertical axis) with the fixed £ = 2, where 

( B ) 

the case of 2Kb"t\s < 1(> 1) denotes the enhancement 
(suppression) of the superfluidity. For N = 0, the 
current-current boson-fermion coupling M suppresses the 
superfluidity. The effect of the suppression reduces and 



vanishes, i.e., 77. 



(B) 



1/(2Kb) on the solid line in Fig. 



[T] due to the compensation effect of these two coupling 
constants. For small \M\ and \N\, r/i B) = 1/{2K B ) is ob- 
tained at M = N for MN > 0, and at M = -N£/(l + £) 
for MN < 0. In the region of N < 0, we find a novel 
feature for small \M\. The superfluidity is enhanced by 
the coupling \N\, and is further enhanced due to the cou- 
pling \M\. The enhancement becomes optimized, i.e., the 

( B) 

minimum value of 2KbT) s is obtained at a certain value 
of \M\. 



sence of fermion-boson interaction only in the limit of 
large fermion density. We have also been able to predict 
how the relative weights of the peaks in Bragg scattering 
intensity depend on the parameters of the phenomeno- 
logical Hamiltonian. 
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Appendix A: Consequences of Galilean invariance 

Let us consider a system with N p different types of 
particles indexed by 1 < a < N p with masses M a . Under 
a Galilean boost, pi >a — > p,^ a + M a v and r^ Q — > r^ a + vt. 
This transformation is realized by the unitary operator: 



U — exp 



N p N, 



2^ 2^ M a vr^ a 



a—1 i a — l 



(Al) 



One has: 



VI. CONCLUSION 

We have analyzed the mixture of bosonic and fcrmionic 
atoms in one dimension using a Green's function equa- 
tion of motion method starting from a phenomenological 
bosonized Hamiltonian. We have derived expressions of 
the zero temperature exponents and of the thermal cor- 
relation thermal correlation lengths for the density wave 
and superfluid order parameters. For the exponents, we 
have recovered the exponents previously obtained for the 
mixture in Ref. [28| and we have studied their behavior 
as a function of the fermion density for fixed interaction. 
W e have found that for weak interaction, the behav- 
ior of the superfluid exponent of the bosons could be 
non-monotonous, although boson superfluid correlations 
were always enhanced with respect to the system without 
fermions. Such behavior can lead to a non-monotonous 
dependence of the superfluid transition temperature as 
a function of fermion density in an array of weakly cou- 
pled mixtures generalizing the bosonic array considered 
in Ref. 49J. This behavior is also in contrast with the one 
observed in experiments with three dimensional interact- 
ing systems [231 ] , where fermion doping was seen to reduce 
superfluidity. In one dimensional systems at incommen- 
surate filling, fermion doping is seen to always enhance 
superfluidity, the maximal enhancement being obtained 
near the collapse instability. By contrast, we have seen 
that density wave exponents are decreasing functions of 
the fermion density, which recover their value in the ab- 



U^HU = H + vP+ 2 My2 > 



(A2) 



where the total mass M = N a M a and the total mo- 
mentum: 



N p N a 



(A3) 



a— 1 i a — 1 



is the total momentum of the system. In second quanti- 
zation, the operator U takes the form: 



„ n p 
i / dry^ y M a vrp a (r) 

"' a=l 



C/ = exp i drX]M a vrp a (r) , (A4) 



and one has : U^il) a (r)U = e lMaVr ip a (r). There- 
fore, in bosonization, a Galilean boost takes the form 
0«(r) -> 0«(r) + M a vr, i.e., 7ril a (r) -> 7rII a (r) + M a v. 
From the conservation equation for particles of type a, 
d x ja + dtp a = and the bosonization relations (IMl) . we 
have that the current density is j a — dt4> a /^- Using the 
equation of motion, we can rewrite j a — ^2 b M a b^8b/'K. 
Therefore, under a Galilean boost, we will have: 



Ua) = -J2 M ab M b v. 

TV < J 



(A5) 



Turning to the particle current, in the rest frame, the 
particle current will be (J a ) = vp a . This can be seen for 
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instance by calculating j a = dt(<f> a (x — vt))/n. From this 
result, it is clear that one must have: 



= J2 MabM b 



(A6) 



for each a. In the case N p = 2 this gives Eq. ([5]). 

An alternative derivation of Eq. (|5|) is obtained by cal- 
culating the ground state energy under twisted bound- 
ary conditions ([3]) using second order perturbation theory 
1501. One finds: 



6E = Y, 



Pa 



2M a L 



-y 

r,2 



(0\J a \n)(n\J b \0) 

E n — Eq 



(p a ip b {A7) 



where J a = P a /M a is the current of particles of type a. 
Because of the translational invariance of the system, the 
eigenstates of the Hamiltonian are also eigenstates of the 
total momentum operator P = ^ a P a . Since the ground 
state is non-degenerate, this implies that (n|P|0) = 0, 
i.e., for the Bose-Fermi mixture: (0|Pjr|n) + (0|PB|n) =0, 



leading to the following expressions for M a b'. 



where: 



M BB 
M FF 

M BF 

X 



E 



M B L M I ' 

kRf 2?r x 

M F L My 

2tt X 



\(n\P B \0}\'< 



(A8) 
(A9) 
(A10) 

(All) 



E n — Eq 

Obviously, x > 0. Using the expressions (|A8j) . we obtain: 
so that det(M) > 0. 
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